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This work presents a novel direct torque and flux control (DTFC) of 
permanent magnet synchronous motor (PMSM) with analytically-tuned 
proportional integral (PI) controllers. The proportional (Kp) and integral (Ki) 
gains of the PI controllers were accurately determined, from first principle, 
using the model of the control system. The PI flux and torque controllers 


were then developed in rotor reference frame. The designed PI controllers, 


together with the torque and flux controllers, were tested on a permanent 
magnet synchronous motor (PMSM). The results obtained were compared 
with results from conventional DTFC system using manually-tuned PI 
controllers. The total harmonic distortion (THD) of motor phase currents is 
18.80% and 4.81% for the conventional and proposed models respectively. 
This confirms a significant reduction in torque ripples. The control system 
PMSM was tested for step torque loading and found to offer excellent performance 
Rotor reference frame both during load changes, speed reversal, and constant load conditions. 
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1. INTRODUCTION 

Permanent magnet synchronous motors (PMSM) have very wide field of application owing to their 
many advantages such as higher efficiency, power factor, power density, and torque-to-inertia ratio compared 
to the traditional induction motor often referred to as the workhorse of the industry [1 ]-[3]. This has resulted 
in intensified research activities in the design, analysis, and control of PMSM. In most industrial applications 
where PMSMs are used, high precision in speed and torque is usually required [4]-[9]. It has been established 
that direct torque control (DTC) offers better dynamic performance compared to field orientation control (FOC) 
especially for spplications requiring high precision, sensitivity, and minimized torque ripples [10]-[14]. 

Advancement of DTC has witnessed much attention in the development of control algorithms 
without considering analytical evaluation of PI controller gains and this has greatly affected precision and 
speed/torque distortion. Different DTC strategies based on hysteresis control were presented in [15]-[17]. In 
these works, inverter voltage vectors were selected from hysteresis-based switching tables but the problem of 
high torque ripples persists. Bo-Wen ef al. in [18], a DTC method based on active disturbance rejection 


Journal homepage: http://ijpeds.iaescore.com 


2104 O ISSN: 2088-8694 


control was presented while [19] presented DTC in flux weakening region. In the hysteresis controllers were 
replaced with PI controllers without any analytical method for selection of the proportional and integral gain 
values of the PI controllers [20]-[22]. The results of these works show various levels of torque ripples which 
are not suitable for applications requiring minimal ripples in torque and speed. Recently, a model-based 
analytical evaluation of PI controller gains for field orientation control (FOC) was developed in [23] and 
tested on a PMSM. The results show that by accurately determining the gains of the PI controller, torque, and 
speed ripples as well as overshoots were greatly reduced. 

The greatest drawback of the traditional DTC is its high torque and flux ripples. However, several 
researchers have tried different methods in an attempt to reduce these ripples and improve the overall 
performance of the traditional DTC. Some of these methods include: the hysteresis regulators-based stator 
voltage vector selection in [16], the active disturbance rejection controller DTC in [18], the MT-frame flux 
weakening DTC in [19], the SVM-based predictive torque control in [20], the PI based DTC with pole 
placement technique in [21], and sliding mode control with PI-based DTC in [22]. Apart from the predictive 
torque algorithm in [21] that significantly reduced the torque ripples, other algorithms did not achieve much 
in terms of torque and flux ripples reduction. 

From the foregoing, it is seen that a clear gap exists in model-based analytical evaluation of PI 
controller gains for DTC systems. In this work therefore, a PI controller tuned by analytical means, using 
mathematical model of the control system, is developed and employed in a novel DTC algorithm. The 
effectiveness of the proposed scheme is demonstrated on a PMSM fed by a space vector pulsewidth 
modulated voltage source inverter (SV-PWM VSI). MATLAB/Simulink was used for modelling and 
simulation in this work. 

This work is organised is being as: section 1 is Introduction while section 2 presents the PI 
controller design developed using the machine qd-axes model. Section 3 contains the design of the PI flux 
and torque controllers while section 4 presents the simulations and results. Section 5 is the conclusion. 


2. PI CONTROLLER DESIGN 
The dq-axis voltage of a permanent magnet synchronous motor can be written in terms of flux 
linkage as [23]. 


dag 


Va = Rig +75- Og (1) 

Vy = Rig + [2+ Ag (2) 
Where 

Aa = Laia + Apm (3) 
And 

eT (4) 
From (3) and (4). 

iq = “em (5) 

u=% © 


Substituting (5) and (6) into (1) and (2) respectively gives. 


Àd- pm da 
Va SiG tae (7) 
— p4a, Wa 
Y= Ri + Git Oda (8) 


Laplace transform of (7) and (8) gives. 
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Va(s) = (È + s)Ag(s) — 044 (8) — “E 


v (8) = (s + [Aq (s) + a(S) 


From (9) 


RApm, 1 


I 


La Ga +s) 


Aa(s) = [Va (s) + @Ag(s) + 


(9) 


(10) 


(11) 


In 11 is represented in block diagram as of Figure 1. It is seen from (1) and (2) that the dq flux control 
loops are not independent due to the back-emf terms in both equations. To make the dq flux control loops 
independent, a back-emf decoupling term is introduced. This decoupling term is introduced just after the PI 
controller to serve as a disturbance with the same value but but having opposite sign as the back-emf term in the 
motor model. The block diagram of the system is now as shown in Figure 2. Since the decoupling term will cancel 
the effect of the back emf term, Figure 2 reduces to Figure 3. There will be both sensor delay and computation time 
delay in the system; a first order time delay is introduced to modify Figure 3 as shown in Figure 4. 


Figure 4. Closed-loop block diagram with time delay 
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The open loop transfer function, Go, is given by (12) as shown in. 


KpSt+Ki 1 La 
s  ‘“Tps+1°Lqs+R 


Go,(s) = (12) 


If the zero of the PI controllers (C K ) is designed to cancel the pole (“Ray La) of the PMSM by pole- 
p 


zero cancellation method, then we let: 


Kip = ko Ia (13) 
From (13), 
K; = K,Kip (14) 


Substituting (14) into (12) gives. 


Kps + KpKip 1 1 


Gor (s) z Š “Tps41 R (15) 
La 
— Kp _ 1 
Gol) = 35 ae (16) 
Let K= a , so that (16) becomes. 
D 
K 
Go, (Ss) = s+ 5 (17) 
The closed-loop transfer function, Gc, is given as. 
PE GoL _ K 
Gc, = 1+Go,  s2+ 7K (18) 
The general equation of a second-order system is [24], [25]. 
oñ 
TUS) = o (19) 
By comparing (18) and (19), we obtain. 
w2 = K or wn = VK (20) 
And 
1 1 
26Wn = Tp or ô = 2TpvK (21) 


Where w,, is the natural frequency and 6 is the damping ratio. The maximum overshoot of a second order 
system is given by [23]. 


Mp = ern (22) 


So, by choosing the maximum percentage overshoot allowed and the value of Tp; Kp and K; can be 
calculated using (20)-(22). 


3. DESIGN OF FLUX AND TORQUE CONTROLLERS 
3.1. PI flux controller model 

In rotor reference frame, the flux of a PMSM is controlled by d-axis current. So, the PI flux 
controller is designed to produce the accurate d-axis reference voltage from d-axis flux error which is its 
input. Thus, to obtain a flux controller model we make use of the d-axis voltage (7). 
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From (7). 


a Aa-Aym daa _ R dia 


: ; da ; ; : 
Let the time rate of change of d-axis flux, aa be equal to d-axis flux error per sample time, a where eg is 
Ss 


the d-axis flux error and T, is the sample time. 


Therefore, 

dag = ea 

wn 24) 
Taking integral of 24 gives 24. 

1 

da = 7 J eadt (25) 

Substituting (24) and (25) into (23) gives. 
= fedt- Apm + -wå 
kma S ea UET; pm + 7 T WAg (26) 
$ Ki _ R 
Since Kj, = ie Ia 
In (25) is re-written as. 
— Ki Ki ed 

Va = ay, J eadt — 7 Apm + 7° — WAg (27) 
In (26) is the required PI flux controller model. 
3.2. PI torque controller model 

The torque equation of a PMSM is given as [23]. 

Te = 1.5pApmig + (La — La )iaig ) (28) 
Where p is the pole pair. 
Since Lg = Lg for non-salient pole machine, then. 

Te = 1.5pApmig (29) 


In (28) shows that electromagnetic torque of a non-salient pole PMSM is directly proportional to the q-axis 
current, but from in (4), 


4q 


Aq = Llai > ig = 7 (30) 
Substituting (29) into (28). 

Te = 1.5pAym—! 1 

e 5 1 9PApm Eg (31) 

Differentiating both sides of (30) with respect to time gives. 

dTe _ 15pApm dåq 

d — lq dt (32) 
Let the time rate of change of torque, z, be equal to the torque error per sample time, e,. So 

Mice (33) 
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Similarly, the PI torque controller is designed to produce the exact q-axis reference voltage from its input, 


which is torque error. Since the torque of a PMSM in rotor reference frame is controlled by q-axis current, 
the PI torque controller model is obtained by modifying the q-axis voltage (in (8)). From (8). 


— pta, 44a 
Vy =H FT + wig (34) 
A 
Let the time rate of change of q-axis flux, “a, be equal to the q-axis flux error per sample time. That is; 


a a (35) 


dt Ts 
åq = - feat (36) 


In (32) and (35) show that the q-axis flux error per sample time is directly proportional to torque error. So, by 


L ; ; eq . A 
4 the q-axis flux error per sample time — is obtained. Therefore, 
1.5pApm Ts 


substituting (35) and (36) into (34), (35) is obtained. 


multiplying torque error by 


Kj e 
Va = ax, J eq At + 7 + wg (37) 


The block diagram of the conventional DTFC as reported in [21] is shown in Figure 5 while the block 
diagram of the proposed DTFC with analytically-tuned PI controller is shown in Figure 6. 


PI Controller 
SVPWM 
Inverter 
PI Controller 


Torque and Flux 


ia 


Observer 


Figure 5. Block diagram of the normal or conventional DTFC system 


Torque and Flux SVPWM 


Controller Inverter 


Torque and Flux 


Observer 


Figure 6. Complete block diagram of the proposed DTFC system 
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4. SIMULATIONS AND RESULTS 

The conventional and the proposed models were modelled and simulated in MATLAB/Simulink 
environment and tested on a PMSM with the parameters shown in Tabel 1. The designed controller 
calculated gains of K, = 326.72 and K; = 31583 while the normal PI controller gains were tuned to K, = 5.75 
and K; = 150 in the simulation. The reference torque varies from 5 Nm to -5 Nm and returns to 5 Nm during 
a total simulation time of 5 seconds. The results show that the proposed model effectively tracked the 
reference torque with almost no ripples while the conventional model produced much torque ripples. 

Figure 7 compares the plot of torque response of the proposed and conventional models. It is seen 
that while the conventional PI controller produced much torque ripples, the proposed model gave almost no 
torque ripples. Another strength of the proposed model is its ability to calculate the gains of the controller 
rather than relying on the error-prone and time-wasting tuning methods used in normal PI controllers. Figure 
8 compares the phase currents of the two models. A difference in phase currents, with that of the proposed 
model being higher and with lower total harmonic distortion (THD) is observed. 

The disparity in phase currents is explained by Figure 9 which shows the i, and ig currents. It is 
seen in Figure 9 that the torque producing current (i,) is essentially the same for both models. However, the 
d-axis current of the proposed model is higher than the normal one. This is because the two models are 
different and the proposed model uses higher gains. So, the control effort pushes any difference resulting 
from the higher gain to the d-axis current in order to ensure that the response effectively tracks the command 
torque. For these reasons, the phase currents which are combinations of d-axis and q-axis currents (by inverse 
Park transform) is different for the two models. The FFT for phase ‘a’ current for the conventional and the 
proposed DTFC are shown in Figures 10 and 11 respectively. The proposed model shows superior 
performance with THD of 4.81% compared to THD of 18.80% for the conventional DTFC. 


Table 1. The models conventional and simulated in MATLAB/Simulink 


Motor parameters Value 
Rated Power 9.4 kW 
Frequency 50 Hz 
Stator resistance (Rs) 0.203 Q 
Constamt rotor flux linkage(Ay) — 0.123Wb 
Inductance d axis (La) 0.0021H 
Inductance q axis (L,) 0.0021H 
Inertia constant (J) 0.0048K gm? 
No. of poles (P) 4 


ia-proposed 
00 f f f F f f ia-normal 
= a 
pan s 
= en l E 
= 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 
time (s) —— ib-proposed 
00r T F l ib-normal 


ib(A) 
Oo 
| 


ic-proposed 
ic-normal 


torque (Nm) 


ic (A) 
= 


0 0.5 1 1.5 2 2.5 3 3.5 4 45 5 
time(s) 


Figure 7. Torque response of the two model Figure 8. Stator phase currents 
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T T T T 
iq-normal : 
mi- ----- iq-proposed |-----;/™ 


< olata a a a aa Coe ee I | Laas nero 
= 
10 1 L L L 1 1 L Ł 1 
o 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 
time (s) 
60 
AO OREN AN: EPESOES i esasa i Tssa ; sed id-proposed |... 
Ea : H í : H id-nomal 
z i : 
o i i i i i i i t i 
o 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 
time (s) 
Figure 9. dq-axes currents 
FFT analysis 
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Figure 10. FFT analysis of phase ‘a’ current for conventional DTFC 
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Figure 11. FFT analysis of phase ‘a’ current for proposed DTFC 


5. CONCLUSION 

This work has presented a novel direct torque and flux control in which the PI controllers were 
tuned analytically from first principle and compared with a conventional DTFC system that employs 
manually-tuned PI controllers. The accuracy of both the control model and the PI controller tuning method is 
evident in the results which show excellent tracking of the reference torque with no overshoot and less 
ripples than the one with normal PI controllers. By being able to calculate its controller gains; the developed 
algorithm saves the control engineer the trouble of tuning PI controllers by trial-and-error method and also 
performs better than conventional PI controllers. The developed model can be adopted for speed or torque 
control of PMSM. The overall objectives of the work have been achieved. 
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